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Abstract: Theories of dark matter that support bound states are an intriguing possi- 
bility for the identity of the missing mass of the Universe. This article proposes a class 
of models of supersymmetric composite dark matter where the interactions with the 
Standard Model communicate supersymmetry breaking to the dark sector. In these 
models supersymmetry breaking can be treated as a perturbation on the spectrum of 
bound states. Using a general formalism, the spectrum with leading supersymmetry 
effects is computed without specifying the details of the binding dynamics. The in- 
teractions of the composite states with the Standard Model are computed and several 
benchmark models are described. General features of non-relativistic supersymmetric 
bound states are emphasized. 
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1. Introduction 



The nature of dark matter is unknown and its relation to the Standard Model (SM) is 
an open question. The recent spate of anomalies in direct detection experiments [1] and 
cosmic ray signatures [2] has motivated re-examining the standard assumptions about 
the identity of dark matter. Most models of dark matter assume that dark matter 
is an elementary particle with no relevant or long range interactions. If supersymme- 
try is present in these models, the supersymmetric mass splittings are so large that 
the supersymmetric structure of dark matter is unimportant. This article provides a 
framework to illustrate the exact opposite case: dark matter is composite with long 
range interactions and supersymmetry breaking effects are small. 

Recent anomalies have several common features that motivate considering dark 
sectors that support bound states. Bound states naturally enjoy a hierarchy of different 
scales. Inelastic Dark Matter explanations of DAMA, e.g. [3-5], require several scales 
to reconcile the anomalies with the null results of other direct detection experiments. A 
hierarchy of scales is also employed in the Exciting Dark Matter scenario [6] to explain 
the 511 keV signal from INTEGRAL/SPI. Additional structure in the dark sector is 
also motivated by positron excesses in cosmic ray data, which might be a result of 
cascade decays in the dark sector. Examining the Standard Model, one finds a variety 
of different bound state systems: mesons and baryons, nuclei, atoms, and molecules. 
Given the prevalence of bound states in Standard Model systems, it is natural to explore 
the possibility [7-9] that dark matter is composed of bound states in a separate sector. 

Fermions with gauge interactions are a ubiquitous ingredient in theories beyond the 
Standard Model. It is plausible that there are additional gauge sectors that SM fermions 
are not charged under. If there are no SM particles directly charged under the new 
gauge interaction, then experimental limits on decoupled gauge sectors are extremely 
weak. If supersymmetry breaking is only weakly mediated to the dark sector, perhaps 
through dark matter's interactions with the Standard Model, then the magnitude of 
supersymmetry breaking effects can be extremely small. This allows for the possibility 
that dark matter is nearly supersymmetric. If there are any bound states in the dark 
sector, the spectrum will exhibit near Bose-Fermi degeneracy. Such weakly coupled 
hidden sectors also naturally sit near the GeV scale, which makes for interesting dark 
matter phenomenology [10] and experimental signatures [11]. 

Investigating nearly supersymmetric bound states arising from perturbative Coulom- 
bic interactions is a relatively intricate process and the standard techniques from quan- 
tum mechanics involve computing first and second order T-matrix elements and then 
diagonalizing the Hamiltonian. At each step the calculation is not supersymmetric, 
although the final answer is supersymmetric. Ultimately, the states have organized 
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themselves into supersymmetric multiplets and the admixtures of different supersym- 
metric particles that each composite state consists of is known. For instance, a spin zero 
fermion-fermion bound state will mix with a spin zero scalar-scalar bound state. Since 
phenomenological applications depend on these admixtures, it would be convenient to 
understand their structure and how they generalize to other bound state systems. Sim- 
ilarly, phenomenological studies would be made easier by understanding how bound 
state interactions are constrained by supersymmetry. This article develops a simple 
formalism to do this using off shell superfields. 

The organization of the paper is as follows. Sec. 2 reviews non-relativistic super- 
symmetric bound states, focusing on how supersymmetry organizes the spectrum and 
superspin wavefunctions of the states. The free effective action is also introduced, 
which will form the basis for computing the supersymmetric interactions of the bound 
states. Sec. 3 incorporates the effects of supersymmetry breaking into the spectrum 
for the case where the dominant source of supersymmetry breaking is the soft masses 
of the scalar constituents. Sec. 4 computes the interactions of the bound states when 
interacting with weakly coupled external gauge interactions. Sec. 5 constructs a re- 
alistic model of nearly supersymmetric atomic dark matter. Sec. 6 discusses possible 
directions of future research for models along these lines, including recombination and 
the formation of supersymmetric molecules. Sec. 7 makes some concluding remarks. 

2. Non-relativistic Supersymmetric Bound States 

This section studies how non-relativistic supersymmetric bound states organize them- 
selves into supermultiplets. Sec. 2.1 outlines a general procedure for determining the 
composition of non-relativistic bound states formed from massive superfields. When 
applicable, this procedure has the advantage of sidestepping a detailed perturbative 
calculation in favor of some superfield algebra. This procedure is illustrated in the 
particular case of bound states formed from two chiral multiplets. Sec. 2.2 continues 
the study of this particular example by introducing an effective field theory descrip- 
tion of the ground state. This will provide the basis for Sec. 4, in which bound state 
interactions are discussed. 

2.1 Wavefunctions From Superspace 

Non-relativistic bound states have a structure that is simple to understand because 
they benefit from a good expansion parameter: the velocity v. This is especially the 
case for two-body systems, where an expansion in powers of v not only helps to organize 
calculations but also determines the relevant scales of the problem. The gross structure 
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of the spectrum can be organized into principle excitations split by energies of order 



m prin oc fiv 2 , 



(2.1) 



where [i is the reduced mass. Fine structure effects are the next order correction in the 
non-relativistic expansion, appearing as 



Recent papers [12, 13] have computed the fine structure of supersymmetric hydrogen 
through explicit calculation. This section rederives these results by considering how 
supersymmetry acts on the bound states in the non-relativistic limit. The organization 
of the spectrum into supermultiplets does not depend on the details of the binding 
dynamics except for specific quantities, such as energy splittings. Most notably, the 
superspin wavefunctions are completely determined by supersymmetry alone if there 
are no accidental degeneracies in the spectrum. This method of using supersymmetry 
to fix the superspin wavefunctions is applicable to a wider class of non-relativistic 
bound states than Coulombic bound states and more cleanly delineates which quantities 
depend upon dynamics versus the structure of supersymmetry. 

For simplicity, assume that the bound state is supported by a central potential 
that is spin-independent at 0(v 2 ). This is true for a wide range of composite states, 
including those bound together by the exchange of light vector or chiral multiplets. 
The ground state then has a non-degenerate radial wavefunction with I = [14] and 
factorizes as 



to leading order, where |«S) is the superspin part of the wavefunction. At leading 
order in v the supercharges act only on \S), leaving \ip(r)) intact, because gradients 
of non-relativistic wavefunctions are suppressed, diip ~ 0(v). Since ip(r) has trivial 
angular dependence, decomposing |«S) into irreducible representations decomposes 
into irreducible supersymmetry representations, ftj, were j refers to the spin of the 
Clifford vacuum (e.g. Qq is the chiral multiplet). At 0{v A ) the Hamiltonian is typically 
spin-dependent and any degeneracy among the Q/s will generically be lifted in the 
absence of any special symmetries. For Q^s that are accidentally degenerate at 0(v 2 ), 
there can be large mixing that depends on the details of the dynamics, though in 
many cases the appropriate mass eigenstates are determined by the action of addtional 
symmetries on the supermultiplets. 

As an example that illustrates this decomposition, consider the model bound state 
system that will form the main subject of this article. It consists of four massive chiral 



7BFS OC /if 4 . 



(2.2) 



|tf > = |V(r)> ® \S) 



(2.3) 



4 



superfields (E, E c , P, P c ) with Dirac masses m e and m p satisfying m e < m p . The 
binding dynamics respect parity, under which the coordinates and superfields transform 

as 

ar" <-> (-l)'V 9 a •H- 9q, P <-> P ct <->■ £ cf (2.4) 

The dynamics also respect a [/(l)R-symmetry and a £7(l) e x U(l) p flavor symmetry. 
The charges of the component fields are taken to be 
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(2.5) 



Significantly, the U(1)r symmetry and Z 2 parity do not commute and combine into an 
0(2)r symmetry. This can be seen by considering the "selectrons" e and e c . Parity, 
V, acts upon the selectrons as 

Ve = e ct (2.6) 

while under a U(1)r transformation, R(a), the selectrons transform as 

R(a)e = e ia e and R(a)e c] = e~ ia e c] (2.7) 

so that [V,R(a)]e ^ 0. Thus U(1)r and Z 2 are not a direct product and instead 
combine as the semi-direct product U(1)r x Z 2 = 0(2)r. This is important because 
0(2)r has two-dimensional irreducible representations that are realized in the bound 
state spectrum. In particular any state that transforms non-trivially under U(1)r must 
sit in an 0{2) R doublet. 

For this system, the superspin wavefunction |<S) in Eq. 2.3 decomposes as two 
chiral multiplets and one vector multiplet (^), as can be verified by counting degrees 
of freedom. As a consequence of the 0(2)r symmetry, however, the two chiral multiplets 
combine into a hypermultiplet Jf? so that the decomposition of \S) reads 

\S) = 2VL ®VLi= ®r (2.8) 

Both Jtf? and Y are charged under the global U(l) e+P flavor symmetry of the theory. 
The superspin wavefunctions of the ground state are fixed by (super) symmetry at 
leading order because J4? and Y are irreducible under the full symmetry group and 
therefore insensitive to mixing. 
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Supersymmetry organizes non-relativistic pairs of free particles into supermulti- 
plets, determining the bound state wavefunctions at leading order in v in terms of the 
constituent particles. The organization of pairs of free particles into supermultiplets 
is found by putting E, E c , P and P c on shell and constructing all possible superfield 
bilinears. The resulting bilinears will have spins ranging from to 1. For example, 
the superfields P and £ ct yield the bilinears PE C \ V a PE c \ PV&E'* and V a PV A E c l 
These bilinears can then be decomposed into irreducible supersymmetry representa- 
tions with the help of projection operators, which in the case of spin zero superfields 
are given by 

V 2f)2 f)2 V 2 VV 2 V 

Vi = V 2 = and V T = (2.9) 

1 16D 16D 8D K J 

where Vi + V 2 + V T = 1 [15]. 

The decomposition is simplified by noting that the same state can appear in many 
different bilinears. In fact the bilinears 

V 2 PE = PE ViP c] E c] = P c] E c] and V T PE C ^ (2.10) 

contain all the states as can be verified by counting degrees of freedom. Expanding the 
first two bilinears in Eq. 2.10 using the non-relativistic fields 1 



gives the superfields (cf. [16]) 

PEcx <f> p (f> e + V2Q a (c^0 e + s (f> P i>°) (2.12) 
-6 2 {s 2 e (t) p (j)l c + (? e 4>pc(j) e - s 2 e(tp p ipe)o S j 

and 

P ct E ct oc ^c0 e c + V2G a (crfitfl + se4<€) (2-13) 

-0 2 (C00 p e c + sl(f)L(j)e + S29(lpp^e)o 



where the dimensionless Q = yjm p + m e 8 a has been introduced, and the mixing angle 
9 is defined by 

iasi 2 e = — (2.i4) 

m p 



1 The superscript D indicates that the spinor is in the Dirac basis, where 70 is diagonal. 
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These two superfields have £/(l)#-charges of ±2 and transform into each other under 
parity; they correspond to the two Qq's in The "V wavefunctions are found by 
decomposing the bilinear VtPE c \ which gives a complex vector (curl) superfield with 
components 



D oc c 2 e(ip P ipe)o + s 2 e{4>lc<p l 

Ai OC Sgtpp4>lc - C e 4>lc1p t 
X 2 OC Se1pp(j)e - Cg^pl/Je 

oc ip p a%lj e 



(2.15) 



Going to the parity eigenbasis and introducing notation for the various states gives 
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C 2 e|(^e)o) + ^|(|0p+0 e -) - |0p-0 e +» 

73(l0p+0e+) - |0p-0 e -)) 

^§(|0p-0e+) - |0p+0e-)) + S20|(^p^e)o) 
= C e \i) p (j) e+ ) + Sel^p+T/'e) 
= Cg\lpp(j) e -) + Se\<j) p -1pe) 
~- 7^(l0p+0e+) + I0p-0e-)) 
= ^g(l0p+0e-) + \<f>p-<Pe+)) 



(2.16) 



where {ip p ip e )o 



), and cg,sg are COS6 1 and sin#, 



respectively. These are the same wavefunctions found in [12,13] by means of a detailed 
perturbative calculation in the particular case of supersymmetric hydrogen. 

Although the Y and Jf? wavefunctions have been determined here without specify- 
ing the binding dynamics, the mass splitting between Y and M 1 can only be determined 
by doing a dynamical calculation. In the absence of any special symmetries, however, 
it is expected that mps = fflf — m^ will be at the fine structure scale, m-Fs ~ 0(v a /j,), 
and in the case of supersymmetric hydrogen one finds 



m FS 



■ 4 



(2.17) 



The states in Eq. 2.16 are organized according to their 0(2) r representations, with 
simple transformation properties under parity, because the breaking of 0(2) r plays an 
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important role in lifting degeneracies in the spectrum once supersymmetry is broken. 
The states \±i £± an d cc± transform in two-dimensional representations of 0(2)r with 
U(1)r charges of 1, 1, and 2, respectively For example, the doublet 




(2.18) 



transforms irreducibly as 

V : zu a 3 zu R(a) : w -»■ e 2iaa2 w (2.19) 

The states v^, and u)± are invariant under R{a) and thus transform as 0(2)r singlets. 

To illustrate the action of supersymmetry on the ground states, consider the heavy 
proton limit, 6 — > 0. In this limit, supersymmetry clocks the states of the heavier 
constituent, leaving the valence particle intact. In particular, the Y states consist of 
a light electron orbiting a heavy proton multiplet and the M' states consist of a light 
selectron orbiting a heavy proton multiplet. 

This method of calculating superspin wavefunctions through decomposing prod- 
ucts of superfields is general and can be applied to a wide class of non-relativistic 
supersymmetric bound state problems. For example, the superspin wavefunctions of 
non-relativistic SU (3) baryons can be found by studying the decomposition of super- 
field trilinears. In this case, acting with the projection operators in Eq. 2.9 on spin 
zero trilinears does not give all of the wavefunctions, and a spin | trilinear is neces- 
sary. Similarly, the study of the bound states of a massive chiral and a massive vector 
superfield requires higher spin projections. 

Excited state wavefunctions 

This prescription for finding the superspin wavefunctions can also be applied to the 
excited states. For a given spatial wavefunction \nl), the various excited states can 
be built by acting with supersymmetry on the Clifford vacua defined by the particle 
content, 

{|nZ)<g> \nl) <g> (a f <g> \Q S )) , \nl) <g> (a j2 \Q s )) } (2.20) 

where \ft s ), a'|fi a ) and a^ 2 1 O s ) are the superspin wavefunctions derived in the previous 
section. For example in the case considered above \Q S ) is either \Q ) or and the 

raising operators fill out the various states in "V and ffl . Decomposing Eq. 2.20 into 
supermultiplets is equivalent to switching to the basis 

{\nl)®\Q s ), a 1 " (8) (\nl) ® , a t2 (\nl) ® \Q S ))} , (2.21) 
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since irreducible representations of super symmetry are obtained by acting with the 
raising operator on Clifford vacua that are irreducible representations of the rotation 
group. This basis switch is just a matter of Clebsch-Gordon algebra and results in the 
decomposition 

Z®n, = fl[j_ s |©...©fl[ l+S |. (2.22) 

For example, in the case considered above, where the bound state is formed from two 
chiral multiplets, the decomposition gives 

/ ® (O g> Qq) = I <g> (fi © tt © Sli ) = fi z _ 1/2 © ^ © ^ © fi l+1/2 (2.23) 

with the two fli related to one another by parity. 

Thus provided that a given Qj does not undergo large mixing, the excited state 
angular/superspin wavefunctions can be found just as for the ground state. One does a 
(single) superfield calculation as before to determine Q s and then transforms from the 
basis of Eq. 2.20 to that in Eq. 2.21 using Clebsch-Gordan coefficients. 

2.2 Effective Action for the Ground State 

Once the ground state spectrum is known, it is important to determine how the various 
states interact with one another as well as with the SM. There are a variety of interac- 
tions, many of which are related through supersymmetric Ward identities. Superfields 
thus offer a convenient method for packaging all these interactions into manifestly su- 
persymmetric forms. This section uses the standard off shell superfield formalism to 
formulate an effective free action for the ground state, postponing until Sec. 4 a dis- 
cussion of ground state interactions. 

J^f is described by two chiral superfields that satisfy the following relations on shell 

Ui oc PE and U\ oc P c] ' E c] '. (2.24) 

A second set of chiral superfields, TL± and T-L^i ls introduced to give the F-terms of %i 
and TV 2 dynamics. The free Lagrangian for M' is given by 

£^ = Jd 4 6 d'^H-H, + n c \ufj + Jd 2 9 ^m H UiU] + h.c. (2.25) 

The equations of motion which follow from Eq. 2.25 then result in the identification 

oc VJ**E and H c 2 oc V 2 P cj E (2.26) 
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With %i and %1 identified as above the appropriate U(1)r and U(l) e+P charges are 
given by 





U 1/ c t 1/c 

rLi 1%2 'M "2 


U(1)r 

U(l)e +P 


2-200 
2 2 2 2 



(2.27) 



so that the Lagrangian is properly invariant under U(1)r and U{l) e+p . Parity acts on 
the composite superfields as 



(2.28) 



so that the Lagrangian is also invariant under parity. 

y is described by an off shell field, V, and an action consistent with the on shell 
constraint V oc VtPE°\ V is a charged vector superfield - a general superfield with 
no Lorentz index 

V(x,9,9) ^V\x,9,9). (2.29) 
The action is written with the help of the supersymmetric field strengths 

W la = ~V 2 V a V and W 2a = ~^V 2 V a V j (2.30) 
which have U(1)r charges of +1. Under parity V and Wi transform as 

V <-> -V and Wf ^ -W 2 a (2.31) 
The free Lagrangian, which is properly invariant under U(1)r and parity, is given by 

C r = Jd 4 9 2m 2 v V j V + Jd 2 9 ^W?W 2a + h.c. (2.32) 

Varying the action yields the equation of motion T> a W" = 2 my V, implying that VtV = 
V on shell. 



3. Supersymmetry Breaking in the Ground State 

The previous section calculated the composition of non-relativistic supersymmetric 
bound states using supersymmetric group theory, focusing on the particular example 
of bound states formed from two chiral superfields. This section builds on Sec. 2 by in- 
corporating the effects of weak supersymmetry breaking on the ground state spectrum. 
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The exact changes to the spectrum resulting from supersymmetry breaking depend 
on the details of the binding dynamics. In many theories, however, supersymmetry 
breaking level splittings induced by the binding dynamics are accompanied by powers 
of the velocity, v. Consequently in the non-relativistic limit supersymmetry breaking 
in the bound state spectrum will be dominated by the differences in the rest energies 
of the constituent fermions and bosons. For such theories the resulting spectrum is 
insensitive to the details of the binding dynamics. 

3.1 Constituent Mass Effects 

The leading supersymmetry breaking effects can be calculated by folding in the per- 
turbed rest energies of the constituents with the ground state superspin wavefunctions 
calculated in Sec. 2. This leading order effect is straightforward to compute if the ef- 
fective scale of supersymmetry breaking in the bound states spectrum, m so ft, is smaller 
than the scale of principle excitations 

m so{t < m prin ~ 0(fiv 2 ). (3.1) 

In this case mixing with excited states is unimportant and the incorporation of su- 
persymmetry breaking into the bound state spectrum reduces to a finite dimensional 
quantum mechanical perturbation theory problem. 

The bound state spectrum has two effective mass scales for supersymmetry breaking 
effects. The first scale is set by the U(l) ^-preserving soft masses, m R _ pres , while the 
second is set by the i7(l ^-violating £?-term masses, uir- Y i \. The breaking of the 
U(1)r symmetry induces splittings between states that are doublets under the 0(2) r 
symmetry. In many implementations of dark sector supersymmetry breaking, [/(l) it- 
violating soft terms will be suppressed relative to the [/(l)^-preserving soft terms and 
for simplicity the relative ordering of the scales is taken to be 

mS-preB^FS > m^-viol (3.2) 

throughout, where mps = 0(fiv 4 ). 

The soft supersymmetry breaking Lagrangian for the chiral-chiral bound state sys- 
tem introduced in Sec. 2 contains a £/(l)R-preserving piece, 

- £ R _ prcs D A e 2 (|e~| 2 + |gf ) + A%\p\ 2 + (3.3) 

and additional supersymmetry breaking terms that break the U(1)r symmetry: 

- £r~ v \o\ ^ B e m e ee c + B p m p pp c + h.c. (3.4) 
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For simplicity the soft parameters are assumed to obey the relations 

A? ~ A? and B p ~ B e = B (3.5) 
In the presence of A~ and B e the selectron mass eigenstates become 

e ± = -L(e±e c t) (3.6) 

with masses 

1 A? 1 

m~ e ± = m e + 5rrie ± = m e + -— ± -5 (3.7) 

2 m c 2 

Analogous expressions hold for the mass eigenstates p±. See Sec. 5.3 for details on a 
particular implementation of supersymmetry breaking in the dark sector that satisfies 
the above assumptions. 

The leading supersymmetry breaking perturbation on the ground state spectrum 
is encapsulated in the perturbing Hamiltonian 

H soit = 5mp±\<j) p ±) (<j) p ±\ + 5m i± \4> e± ){(j) e± \ (3.8) 

The [/(l)R-preserving contributions of if so ft will appear in the combination 

m so ft = ^(Srrip + + dm^ + 5m i+ + <5m e -_) (3.9) 

The rest energy perturbations can now be read off directly from the supersymmetric 
wavefunctions in Eq. 2.16. For example consider the state 

k-) = c 2 e\('4'pipe)o) + S2e(|0 P +0e-) - l0p-0 e +))/v / 2 (3.10) 

The fermion-fermion component is insensitive to H so f t , but the scalar-scalar component 
results in a perturbation 

Am f _ ~ (<j_|# soft k_) = -s\ {5mp + + 5m- e _ + 5m^ + 5m e+ ) = s 2 2e m soft (3.11) 

which is the leading supersymmetry breaking contribution to the mass of in the 
limit that m FS ^> m soft . For many physical applications, such as decays or scattering 
off of SM nuclei, knowing only the leading breaking is sufficient. Using the superspace 
approach for finding the wavefunctions, as in Sec. 2.1, the leading supersymmetry 
breaking can thus be found for a broad range of perturbative bound states. 
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3.2 Subdominant Effects 

Supersymmetry breaking effects begin to grow in complexity beyond the rest mass 
perturbation. The next most important term in the non-relativistic expansion is the 
kinetic energy perturbation 

- -t- < 3 - 12 > 

This changes the principle structure of the bound state and leads to a 0(v 2 ) perturbing 
Hamiltonian 

H v 2 = - (^-^j (cos 4 # 5m g ±\(f) e ±)((f) e ± \ + sin 4 e 5m p ±\(f) p ±) (<j) p ±\) . (3.13) 

At the level of fine structure many new effects arise. These include additional kinematic 
effects from 0(p 4 ) terms and, in the case of supersymmetric hydrogen, gaugino mass 
effects and D-term contributions. Incorporating all these effects requires using the T 
matrix and computing all tree-level Feynman diagrams contributing to ep — > ep matrix 
elements. The T matrix is proportional to an effective non-relativistic Hamiltonian 
that can be used to do perturbation theory, as in the calculation of the fine structure 
of supersymmetric hydrogen [12,13,17]. 

3.3 Eigenstates 

In this section the ground state spectrum with weakly broken supersymmetry is pre- 
sented by diagonalizing the perturbation H so{t . 

Scalars 

In the absence of supersymmetry breaking the hypermultiplet contains the degenerate 
pair of positive parity scalar bound states zu + and u> + . In the presence of H so{t these 
states mix maximally: 



ro + V ( m soh B \ 



I "'■soft J-> \ I zu + 

V B m so{t ) \u + 



(3.14) 



B characterizes the size of 0(2) r breaking and mixes states of different i?-charge. The 
mass eigenstates are 

ou 1+ = 0p+0 e + = —t=(^+ + w+) m Wl+ = 5me+ + 5mp + = m so{t + B 

V2 

^2+ = P -0e- = —j={^+ - w+) m W2+ = Srrie- + 5m,~ p _ = m soft - B. (3.15) 

a/2 
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2 m soft S 29 
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2 m soft S 29 
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2 TO soft S 29 
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Figure 1: The ground state spectrum for the case ttt-fs > and in the limit that mps S> w so ft, 
where the mixing between the vector multiplet and hypermultiplet is small. The composition 
of the various states (see Eq. 2.16 and Eq. 3.15) is a function of tan 2 9 = m e /m p . For 
example, in the heavy proton limit (m e <C m p ), U)- is predominantly scalar-scalar, whereas 
in the opposite limit (m e ~ m p ), cj_ is predominantly fermion-fermion. 



In the supersymmetric limit the ground state contains three parity odd scalars, 
one of which, <j_, is in the vector multiplet and two of which, io_ and cc7_, are in the 
hypermultiplet. In the presence of -ff so ft and H v i all three states mix: 

c_ \ / m FS + s| m so ft -|s 4e m soft \s Ae B' \ / \ 
w_ I -^ 4 em S oft C2 m soft -c 2 ,^' wJ (3.16) 

Here B' characterizes the £/"(l)fl-breaking in this sector and comes about through H v % 
in Eq. 3.13 or from the difference in the 5-term masses between the e and p and is of 
the order 

B' ~ 0(Bv 2 ), 0(B e — B p ). (3.17) 
By specializing to the regime where 

B' < m FS , tan 2 26>m soft (3.18) 
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so that mixing between the three states occurs primarily between and u_, one 
obtains simple formulae for the approximate energy levels: 



«1fs + ra S oft / m FS - c 4e m soft \ / s\ e m 2 solt 

2 V 2 /V ( mp s _ c 4 0m soft ) 2 



m ^ = | ""sort + ^-Ban: M ^ ^sott (3^ 



"'i-s +- '",.,1. / m FS - c 4e m soit \ I | 4? m soft j;5 2()) 



2 V 2 J V ^ mFS _ C 49 m soft) 2 

= "Vft (3.21) 

Non-zero 5-terms split the 0(2)r doublet containing w + and vj— 

In the limit m so f t /mFs — > oo with < 6 < | (respectively | < 9 < |), the state 
(respectively w_) becomes the (i/j e i/j p )o bound state. Naively, for 9 2 ~ ^ ~ and 
m FS — | a vA t the splitting between this state and the vector (ip e tp p ) 1 should give the 
hyperfine splitting in regular hydrogen; however, Eq. 3.19 yields instead: 

2 2atml 

m (ip e 4> P h ~ m (^ P h = m v» ~ "V ->■ s 2g m FS -> (3.22) 

Trip 

This is not the correct hyperfine splitting of regular hydrogen which is 

8 a 4 ml 

muFS = 7j 3.23 

3 m p 

for a point-like proton. This difference arises because m so f t /mFs — )■ oo is not the full 
decoupling limit. In particular, the ground state of supersymmetric hydrogen contains 
admixtures of higher principle excitations arising from gaugino exchange at second 
order in perturbation theory. These effects contribute to the hyperfine splitting in 
Eq. 3.22 but disappear in the full decoupling limit where the gaugino mass goes to 
infinity, my — > oo. 

Fermions 

In the absence of supersymmetry breaking the vector multiplet (hypermultiplet) con- 
tains the degenerate pair of j — | bound states £± (x±)- I n the presence of H so f t the 
states of equal parity mix with one another: 

T 

fX±\ f m FS + a S 20 m soft ±2 B -\s 4e m soit \ fx±\ /g 

V 6fc / V -\s Ae m soit (4 + 4)m soit ±\BJ \£± ) 
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The spectrum is given by 

m FS + m soft / m FS - 4 e m so{t \ 



m 



x± 



2 + V 2 J 

m FS + m soft / m FS - 4 9 m so{t \ 




12 2 
4 S 40 m soft 



±-£ (3.25) 



(m FS - c§„m soft ) 2 2 



For non-zero B the 0(2)r symmetry that ensured the degeneracy of the pair of states \± 
as well as the pair of states £± is broken and the fermionic spectrum splits completely. 

Vector 

The vector state \v) = |(Vy0e)i) is insensitive to H so ^ and, as a consequence, does not 
feel supersymmetry breaking at leading order. 



4. Interactions 

Composite systems have a wide range of interactions that are controlled by selection 
rules and form factors that results in these systems having a much richer phenomenology 
than elementary particles. This section uses the effective field theory of Sec. 2.2 to study 
the interactions that arise when composite states inherit gauge interactions from their 
constituents (cf. [18]). 

Sec. 4.1 considers the case where the constituents are charged under an unbroken 
vectorial gauge symmetry U{l)y such that the composite state is neutral with the 
following charge and parity assignments 





E E c P P c 


U{l)y 


-1 +1 +1 -1 



and 



Vh-V 



(4.1) 



V does not need to be responsible for binding the chiral multiplets together; e.g., 
the binding could arise from a Yukawa force. The U(l)v gauge interactions of the 
constituents induce a number of effective operators, including charge radius, Rayleigh 
scattering, and spin flip operators. Specializing to the case where the hypermultiplet 
H is lighter than the vector multiplet V and m F g 3> m so ft, decays within the ground 
state are discussed in detail. It is found that the states of the vector multiplet V decay 
relatively quickly down to H, while the decays within Ti are much slower. 

Sec. 4.2 briefly considers the case where the constituents are charged under a broken 
axial gauge symmetry £/(1)a with charge and parity assignments 





E E c P P c 


U(1)a 


+ 1 +1 -1 -1 



and 



A A 



(4.2) 
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In models such as that of Sec. 5 where A undergoes kinetic mixing with the SSM, these 
interactions mediate the dominant coupling of dark atoms to the Standard Model. 
Sec. 4.2 discusses the allowed scattering channels and finds the leading supersymmetric 
axial interactions. 

4.1 U(l)y interactions 

The interactions of neutral bound states with an external vector superfield, V, are 
characterized by two scales corresponding to the charge radius, R e , and magnetic radius, 
R m . Physically R e corresponds to the size of the bound state, R e ~ y (r 2 ) . In the case 
of Coulombic bound states R e is given by the Bohr radius, R' 1 = ayy. R m is just the 
Compton wavelength, i?" 1 = /i. For convenience, this section will restrict its discussion 
to supersymmetric hydrogen, although it is generally applicable to chiral-chiral bound 
states. 

Before considering the supersymmetric case, it is instructive to review the leading 
interactions of the photon with the spin-singlet ground state of regular hydrogen. The 
leading elastic interaction comes from the charge radius operator 

SVCM^(W e )^(W e )o^ (4.3) 

which is fully determined by the charge distribution of the bound state. The leading 
inelastic interaction comes from the magnetic spin-flip operator which is determined 
by the fermion content: 

gvRrnd^MiA^eW^ (4.4) 
Finally there is the Rayleigh scattering operator 

g v Rlm H (^ e )l(^ e ) V^ (4.5) 

which makes the sky blue. All other operators are higher order in either gy or yT 1 . 

The next step is to find the set of operators necessary to satisfy the supersymmetric 
Ward identities. V <H> V interactions can be important for scattering processes if the 
states of V are long-lived but have a subdominant effect on the lifetimes of the states 
in V. Because the leading V <H- T-L decay is relatively fast V tends to be short-lived and 
therefore V <H- V interactions are ignored here. The H ■H- H interactions are found in 
the following, since they determine the relaxation timescale of the H supermultiplet. 

The charge radius operator in Eq. 4.3 only depends on the charge distribution, and 
therefore the scalar-scalar bound states must share identical (diagonal) interactions: 



5f V C20 J Re(0p±0e±) t ^(0p±0 e ±)^V 



(4.6) 



Rewriting Eq. 4.3 and 4.6 in terms of the wavef unctions in Eq. 2.16, the charge radius 
interactions become 

g v c 2e R 2 e (4<V^± + ^±d^ ± + dM^ (4.7) 

Similarly, the spin-flip and Rayleigh scattering operators become 

gyRmdM^- + sjjw-)^ (4.8) 

and 

gvR 3 e (m e + m p ) U±uj± + w±w± + sU^V^V*", (4.9) 

respectively. The operators in Eq. 4.7 to 4.9 represent the leading single photon and 
two-photon interactions for the scalar states in V and %. Several interactions remain 
to be found, e.g., the leading single photino interactions as well as the interactions 
for the fermionic states. The coefficients of the remaining interactions are found by 
forming operators from the effective fields of Sec. 2.2. The matching coefficients are 
determined by expanding the supersymmetric operators in terms of their components 
and identifying the corresponding interactions from Eq. 4.7 to 4.9. This procedure 
allows for the various supersymmetric interactions to be systematically enumerated by 
building upon the known interactions of regular hydrogen. 

Interactions of the Hypermultiplet with Higher States 

A variety of processes cause the decay of the excited states to the ground state. For 
example, supersymmetric hydrogen inherits the (fast) electric dipole and magnetic 
dipole transitions of regular hydrogen. Decays from V to "H, however, are not as fast 
and merit further discussion. 

The states in V are connected to % through two one-photon operators of dimension 

five: 

C vw = CMS 20 gvR m Jd 2 e(n c 1 W 1 + U c 2 VV 2 )Wv+ h.c. (4.10) 

+c' M (gv, 6)gvR m Jd 4 6 {U\ + H c 2 )VVWv + h.c. (4.11) 

These two operators are the most general forms for H, V interactions mediated by 
£7(l)v- Higher dimensional operators can be reduced to these two forms with additional 
factors of d 2 acting on Wv by using the matter field equations of motion. 

Only Eq. 4.10 contains the magnetic spin-flip interaction, uj-d ll v u Fy V . The factor 
of s 2 q is fixed by comparison with Eq. 4.4. In supersymmetric hydrogen, some of the 
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component interactions contained in Eq. 4.10 arise from 0{a.y) mixing between the 
ground state and higher principle excitations. For example, excited u + states (2p, 3p, 
etc.) mix with allowing for w M to decay to u + through electric dipole transitions. 
This mixing with excited states is the origin of the "electric" interaction uj + d ll v v Fy 1 ' 
contained in the operator of Eq. 4.10. In this sense, the operator of Eq. 4.10 is nei- 
ther purely magnetic or electric. The role that excited state mixing plays in ensuring 
this supersymmetric result is familiar from the calculation of the supersymmetric spec- 
trum in [12,13,17], where second order perturbation theory is needed to determine the 
spectrum to 0(a v ). 

The operator in Eq. 4.11 does not mediate decays in the supersymmetric limit. 
This can be seen by using the equations of motion to replace VWy with the current Jy. 
Decays through this operator are kinematically forbidden because the mass splitting 
between V and % is much smaller than the mass of any particle charged under U{l)y. 
For this reason we leave the coefficient of this operator undetermined, noting however 
that it can only come in at higher order than g\jR e , since it contains off-diagonal scalar- 
scalar transitions, which do not arise from charge radius scattering. 

The various decay channels induced by the interactions in Eq. 4.10 cause each state 
in V to have the same inclusive decay width to the states of % in the supersymmetric 
limit — otherwise the component propagators of V would have different poles. Therefore, 
the decay width can be calculated by considering the state with the simplest decay 
modes, in this case 

C vm/ D c M g\/FLmm\] ^= (iff 75 + £-) Av <J- (4.12) 
2v 2 

Here Ay is the four-component Majorana gaugino of f (l)v an d rn v is its mass. This 
gives the decay rate 

2 

IVkhv ^ \c M gvS29R m mv\ 2 , F = \c M \ 2 a^ii (4.13) 

Anrriy 

This is a factor of faster than the corresponding spin-flip transition in regular 
hydrogen, which scales as am| s . This is because the decays are dominated by A v 
emission rather than V M emission, for which the amplitude carries an additional factor 
factor of E*, where E is the energy of the emitted gauge particle. 

Interactions within the Hypermultiplet 

Supersymmetry restricts the form of possible interactions significantly, and these re- 
strictions are particularly severe for interactions connecting two chiral superfields. For 
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instance, the only allowed single photon operator, up to possible additional factors of 
d 2 , is 

J d A 6 $i$jj£>>Vv + h.c. (4.14) 

In the case of H •<->■ H interactions, the only operators of this form allowed by the 0(2)r 
and U(l) c+P symmetries of the theory are 

Jd*0(HiH\ - n 2 nl)VWv and J <Pe{UlU^ -UlU^VWy (4.15) 

These operators contain terms like w' i + d fl w + d u Fy 1 ' and u}+d fJi u + d u Fy U , respectively, and 
thus correspond to charge radius interactions. Matching to Eq. 4.7 then gives the 
supersymmetric completion of the charge radius interactions: 

Amv = c E g y c 2e R 2 e ^dH{U\U x - U\U 2 - U\W$ + H c 2 ] H c 2 )VWy (4.16) 

Replacing DWv with the current J\j gives atom-ion scattering. Similarly matching 
onto the Rayleigh scattering operator in Eq. 4.9 yields 

Chhw = c' E gyRl J d 4 6 {UxU\ + ^K^WvWv + h.c. (4.17) 

Just like the operator in Eq. 4.16, this operator will mediate decays within the hyper- 
multiplet once supersymmetry is broken. 

The restriction to operators of the form in Eq. 4.11 is a supersymmetric analog of 
the statement that any interaction involving two scalars and a field strength can be 
written as "(derivatives) x ipd^cj)' d v F^ u ," which implies that transitions between scalar 
states cannot proceed via single photon emission. Thus, for example, direct single 
photon/photino decays from the 2s hypermultiplet to the ground state hypermultiplet 
are forbidden. The decay will instead proceed through either two photon/photino 
transitions or a cascade decay via magnetic operators of the form in Eq. 4.10. 

Hypermultiplet Decays 

In the supersymmetric limit, the hypermultiplet is exactly stable. Once supersymmetry 
is broken and decays within the hypermultiplet become kinematically allowed, it is in- 
teresting to ask what decay channels determine the relaxation timescale. This question 
is complicated by the fact that supersymmetry breaking enters the physics of decays in 
a number of ways. On the one hand, supersymmetry breaking perturbs eigenvalues and 
eigenstates; this opens up phase space, changes the equations of motion, and induces 
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decay channels through mixing. On the other hand, supersymmetry breaking perturbs 
the effective interactions of the non-relativistic constituents. The rest of this section 
considers these possibilities in more detail, with the conclusion that eigenstate mixing 
in the magnetic spin-flip operator, Eq. 4.10, induces the largest decay rates. 

In the presence of soft masses, the supersymmetric operators in Eq. 4.16 and 4.17 
can mediate decays within the hypermultiplet. In the case of the three-body decays 
mediated by the Rayleigh scattering operator in Eq. 4.17, these soft masses appear in 
the eight powers of phase space: 



e 

27 



647r 3 m% 



„' 1 2 



1 fil 16tt 2 V 2 J \m FS J 

In the case of the two-body decays mediated by the charge radius operator in Eq. 4.16, 
these soft masses appear in phase space as well as in an overall factor of m?. This 
latter factor arises from the modified equations of motion for Ay, which imply that 
VWy 3 9$ Ay ex m v #A v . The resulting decay rate is 

r{ Wl+ , W2+ , ro _}_> ?± A v ~ r^-^Av ^ \c E g\jc 2 eRl\ 2 \my(m H -m ?± )| 2 

- \ce \c 2e — [ — ) — ) fi. 4.19 

16 V 2 / V m FS / \ m FS/ 

Here two powers of m so ft arise from cancellations between the terms involving % and 
"H c in Eq. 4.16. Higher order operators may not have this cancellation. 

Next consider how the magnetic spin-flip operator in Eq. 4.10 induces decays in 
the presence of supersymmetry breaking. Mixing between, e.g., the fermionic states x± 
and £± allows all the states in Ji to decay down to through Eq. 4.10, which contains 
interactions of the form 

f 



c M gvR m S29myA v \^JX (J^j (4-20) 



Comparison with Eq. 3.24 shows that, in Eq. 4.20, this fermionic mixing is accounted 
for by making a replacement of the form 

s Ae m soit . ^ 

X^X + ~ A £ 4.21 

4 m FS 
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which leads to the decay rate 



\CMgyS2eR m my\ — 

V 4 m FS / 



> -^g^-soft) 



Anxra 



n 



, 129/ s 2eS46»\ 2 
\Cm\ Oty ' 



(^)\ (4.22) 

V^FS / 



V 8 J 

This decay rate also receives contributions from supersymmetry breaking in the effec- 
tive Yukawa operators of the non-relativistic theory, since the coefficients carry factors 
of m e and m p from the non-relativistic normalization of the scalar constituents 
in Eq. 2.11. These contributions, however, are parametrically smaller by an amount 
0(mp S /ml). Hence the decay rate Eq. 4.22, which is suppressed by only four pow- 
ers of the largest supersymmetry breaking spurion, m so f t , characterizes the relaxation 
timescale of the hypermultiplet. 

4.2 U(1) A interactions 

This section outlines the dominant interactions between dark atoms and an axial U(l) 
with charges given in Eq. 4.2 and mediated by a vector superfield A. Axial gauge 
symmetry forbids mass terms for fermions and therefore the gauge symmetry must 
be broken if non-relativistic bound states exist. As in the previous section, there are 
several allowed supersymmetric operators and the interactions of the vector boson are 
sufficient to fix the coefficients of the operators. 

The leading A M interactions are determined by the axial charges of the constituents. 
The scalars w± in Eq. 2.16 have zero axial charge, but the combinations 

(u + =F c 2 0W_ ± (4.23) 
have charges of ±2 respectively. This leads to the inelastic interactions 

2ig A (c 26 u\Xu~ ~ sW+tys-) A M . (4.24) 
Similarly, the fermion-fermion bound states are charged with interactions given by 

g A m v (c 26 d + s 2 ewl) u^A". (4.25) 

The interactions of Eq. 4.23-4.25 can be embedded in the following superspace opera- 
tors: 

C nnA oc g A c 29 Jd A 9 (U C }U\ + U$u£) A (4.26) 
C WA oc g A s 26 my [d 4 6 {U\ - n 2 r )VA + h.c. (4.27) 



£vva oc g A mlc 2e Jd 4 6 V f VA. (4.28) 
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The interactions of the linear superfield eaten by A can be obtained by going out of 
unitary gauge 



In models like that of Sec. 5, where A M undergoes kinetic mixing with Standard 
Model hypercharge, A M interactions mediate the dominant coupling of dark atoms to 
the Standard Model. This setup can be used to realize inelastic dark matter because 
the elastic interaction of the ground state, U- , with A M is forbidden due to parity. After 
supersymmetry breaking, ^(^+ ± become mass eigenstates. The interactions in 
Eq. 4.24-4.25 then allow w_ to upscatter to ui + and 002+, which are heavier by an 
amount ~ 0(m so { t ), and to v^, which is heavier by an amount ~ 0(rriFs)- 

Higher dimension operators also contribute to the interactions with the standard 
model. For example, the axial spin flip operator 



leads to scattering which can be important in certain regions of parameter space. 

Two body decays mediated by A M are either kinematically forbidden or severely 
suppressed, since m A can only be made smaller than m FS by choosing g A < O(oy). 
Similarly, three body decays mediated by an off-shell A M are subdominant. 

5. Kinetically Mixed Supersymmetric Hydrogen 

This section constructs a minimal model for a nearly supersymmetric dark sector that 
supports Coulombic bound states. Sec. 5.1 introduces a minimal Higgs sector and 
discusses how kinetic mixing of the dark £7(1)a with hypercharge in the supersymmetric 
Standard Model (SSM) drives gauge symmetry breaking in the hidden sector. Sec. 5.2 
adds matter fields that are charged under a second Abelian gauge symmetry, U(l)y, 
that introduces hydrogen-like bound states into the spectrum of the theory. In the 
low energy limit this theory reduces to supersymmetric QED with two massive flavors. 
Sec. 5.3 discusses how supersymmetry breaking is communicated to the dark sector from 
the SSM. Finally, Sec. 5.4 illustrates the scales of the resulting model by calculating 
three benchmark points. 

5.1 Kinetic Mixing 

Abelian field strengths are gauge invariant and therefore no symmetry principle forbids 
mixed field strength terms [19]. Kinetic mixing occurs in extensions of the Standard 



A^ A + 



TTa + TTa 

\/2mA 



(4.29) 




(4.30) 
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Model with additional U(l) gauge factors if there are fields that are charged under 
both the new U(l) and hypercharge. In supersymmetric theories, the entire gauge 
supermultiplet undergoes gauge kinetic mixing, leading to both gaugino kinetic mixing 
and D-term mixing [21,22]. If there are light fields charged under the new U(l), then 
kinetic mixing drives gauge symmetry breaking. 

Consider a minimal example where a dark U(1) A couples to a pair of chiral su- 
perfields $ and $ c with charges ±2 (chosen for later convenience). The Lagrangian is 
given by 

^Hidden = Jd A 6 ($t e 4 ^ A $ + $ c t e - 4 9AA $C ) + JtfQ (lyyjJ _ + Wq) + hx . ( 5 .1) 

where A is the supersymmetric gauge potential of the hidden U{l) Al Wa is the su- 
persymmetric gauge field strength of A, and Wy is the supersymmetric gauge field 
strength of SSM hypercharge. For e C 1 the hidden sector is only a small perturbation 
to the SSM so that all SSM fields will have their normal vacuum expectation values; 
in particular the SSM Higgs fields will acquire vevs along a non D-flat direction: 

2 

D Y = cos 2/3 (5.2) 

This SSM vev now acts as a source term for Wa in Eq. 5.1 and forces <j), the lowest 
component of $, to acquire a vev, since Dy acts an effective Fayet-Illiopoulos term for 
U(1) A . The resulting effective Lagrangian is 

C D = -lDi + D A (eD Y -2g A (\<j ) \ 2 -\ ( pr)) => |0| 2 = |0f + ^ ^ (5.3) 

This D-term potential has a residual flat direction, which can be lifted by Wq. This 
section uses a superpotential 

W = A5($$ c - nl) (5.4) 

where S is a new singlet chiral superfield. With the addition of Wq the vevs of all fields 
are fixed and there are no massless fermions. It is convenient to let the superfields 
acquire vevs and to expand around the new field origin 

($) = v A cos(3 A ($ c > = v A sm(3 A ($<K C ) = $ (5.5) 

where the last expression is enforced by the F-term for S. Solving for tv\ and tan /3a 
gives 
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In the limit /Xq <C eD Y /g A , tan/3 A — > and in the opposite limit, tan/3 A — > 1. Fluctua- 
tions around the vacuum in Eq. 5.5 can be diagonalized using the field definitions 

$ = ( Va + tt) cos /3 a + sin f3 A S c $ c = (v A - vr) sin f3 A + cos /3 A 5' C (5.7) 

so that the superpotential becomes 

W = AS($$ C - nl) = Xv A SS c + ■■■ (5.8) 

clearly showing that S picks up a Dirac mass mj = Av A . In the Kahler term, the 
super-Higgs mechanism takes place: 



K = $V9aA $ + $c t e -4 gA A $c + 

= + S^S + m A (^A + I+J^.) 



2 



+ ••• (5.9) 



The superfield it is clearly identified as the eaten linear superfield, and the vector field 
has picked up a mass 



m A 



2y/2g A v A . (5.10) 



In addition to driving U(1) A gauge symmetry breaking, kinetic mixing also leads 
to £/(l)i{-breaking mass effects in the gaugino sector of the theory. This is because for 
e 7^ gaugino kinetic mixing between U (l)y and U (1) A entangles the U (l)^-preserving 
Dirac mass m A with the f/(l)^-breaking bino mass M\. The effective Lagrangian for 
the gauginos is 

C x = Xyifi^Y + A A z^A A + xJ^Xtt - (MiAyAy + e\ A 0\ Y + m A x^X A + h.c. ) (5.11) 

where Xn is the fermion component of the linear superfield it. The three eigenvalues to 
0(e 2 ) are 



m = mA (l-^0, - m J^^-\, M.fl-^KW) 
Mi — m A J \ Mi + m A J \ mf — m A ' 



Notice that the two mass eigenvalues at \m\ ~ m A are no longer identically the same 
due to kinetic mixing with Ay, and this introduces U(1)r breaking into the hidden 
sector. 

5.2 Charged Matter 

This section adds light, charged matter to the dark sector. The charged matter consists 
of four chiral superfields E, E c , P, and P c that have axial charges under U(1) A . The 
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charge assignments of the dark electron and proton are chosen to be chiral to prevent 
them from acquiring supersymmetric masses in the absence of gauge symmetry break- 
ing. Once the {7(1)a gauge symmetry is broken, these states acquire masses at a scale 
set by mA and tan /3a- In addition to new matter fields, the gauge sector is extended 
by a second gauge group U (l)v under which E, E c , P, and P c have vector-like charges 
and which will lead to the formation of hydrogen-like bound states in the hidden sector. 
In summary the additional matter content has the following charge assignments: 





E E c P P c 


U(1)a 


-1 +1 +1 -1 
+1 +1 -1 -1 



(5.13) 



The superpotential in Eq. 5.4 is augmented by Yukawa terms 

W = W + WV ukawa W Yukawa = y e $ c EE c + y p §PP c + h.c. (5.14) 
so that after £7(1) a breaking both E and P acquire Dirac masses 

WW. = m c (l + Zl±^M! ) EP + m ^ 1 + !L±i^M!) PP' + h.c. (5.15) 
where 

m e = yeSmfikVb and m p — y p cos /Sa^a (5.16) 

None of the fields charged under £7(1) v acquires a vev, and therefore U(l)y is a massless 
gauge multiplet. 

The interactions of the U(1)a vector superfield A with the matter superfields are 
given by 

K = 2g A (a + !LtM (£t£ + E c\ E c _ p\ p _ P c\pcy ( 517 ) 

Here the interactions of the ir have been moved from the superpotential to the Kahler 
potential with the equations of motion. The ir fields can have subdominant mixing 
with the Higgs fields of the SSM and mediate subdominant interactions. 

The Higgs trilinear coupling, A, in Eq. 5.4 and the axial gauge coupling, g A , are 
taken to be large enough that the masses in the axial/Higgs sector are of order aym e or 
larger. With this choice of parameters the axial and Higgs sectors decouple, and the low 
energy limit of the theory is supersymmetric QED with two massive flavors and weakly 
broken supersymmetry. Because the axial/Higgs sectors respect the 0(2)f> symmetry, 
the arguments of Sec. 2.1 go through and, in particular, the leading order superspin 
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wavefunctions are as given in Sec. 2.16. The dominant residual effect of the axial/Higgs 
sector is to perturb the mass splitting between the hypermultiplet and vector multiplet. 
These contributions are suppressed through a combination of coupling constants and/or 
Yukawa suppression. Although the axial U(1)a gauge sector plays a subdominant role 
in the internal dynamics of the hidden sector, it mediates the dominant coupling to 
the Standard Model. In particular it mediates supersymmetry breaking, which is the 
subject of the next section. 

5.3 Supersymmetry Breaking 

Although the hidden sector is supersymmetric at tree level, at the loop level small 
supersymmetry breaking effects are induced through the kinetic mixing portal to the 
SSM. This section discusses the strength with which the constituent particles' masses 
feel supersymmetry breaking. These soft masses determine the leading supersymmetry 
breaking effects in the ground state spectrum, as discussed in Sec. 3. 

The soft parameters to be calculated (see Eq. 3.3 and Eq. 3.4) are the U(1)r- 
preserving Af and A| and the U(l) ^-breaking B e and B p . The largest soft parameters 
are the £7(l)#-preserving ones. If supersymmetry breaking is mediated to the SSM 
through gauge mediation, then these are given by 

A^A^^M? C (5.18) 

where M^ c is the SSM right handed selectron mass. The next largest soft parameters are 
the [/(l)^-breaking I? M -type terms, B e and B p . To isolate how [/(l)^-breaking effects 
are mediated from the SSM, it is useful to integrate out the bino, which generates the 
operator 



e 2 M 1 U 
' □ + M\ 



O t = A A m ; 1 A A + h.c. (5.19) 



B e and B p are then generated upon insertion of 0$ in a loop, with a logarithmically 
enhanced contribution that is the same for both, 

B e ^B p ^B = -A— -log-^ (5.20) 

where Ajjy is the messenger scale. The 5-terms feed into the U(l)y gaugino mass, 
which is highly suppressed due to the indirect communication of £/(l)R-breaking 

a y B a A a v e 2 
m v ~ — A — ~ i a M i- (5.21) 

V 47T (47f) 2 
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In this model, Ay is always light and its mass is smaller than the level splittings induced 
by supersymmetry breaking, which are of order m so ft 



^ ~ M|7 « (5 - 22) 

This justifies ignoring the contributions from niy to the ground state energy levels. 

Sec. 3 described how the dominant communication of supersymmetry breaking to 
the spectrum is through supersymmetry violating perturbations to the rest energies of 
the constituents. Supersymmetry breaking also introduces several dynamical contribu- 
tions to bound state spectroscopy from the exchange of particles frfom the axial U (1)a 
and Higgs sectors. However, these contributions are suppressed for the same reason 
that the supersymmetric contributions from the axial and Higgs sectors are suppressed. 

5.4 Benchmark Models 

This section constructs three benchmark models to illustrate the scales that emerge in 
the hidden sector. Although doing detailed direct detection phenomenology is outside 
the scope of this paper, in both cases we aim to construct spectra compatible with iDM 
phenomenology. In particular we require that the bound states have a mass m^M ~ 
100 GeV and a splitting 5 ~ 100 keV between the ground state and the next highest 
state accessible through axial scattering. Furthermore the (predominantly inelastic) 
scattering cross section between the ground state and standard model nucleons should 
be C(10- 40 ) cm 2 . 

It is possible to meet these criteria; however, some tension exists between meeting 
all three criteria simultaneously. In models consistent with these requirements u;_ or 
is the lightest state. The three states to which can upscatter by exchanging 

an axial photon with a nucleon are u)i+, 0J2+, and v^. The cross section for w_/o"_ to 
upscatter to is velocity suppressed and can be ignored in the following. The cross 
section for cu- to upscatter to u\ + or U2+ via the operator in Eq. 4.24 is given [8,22] by 

RA «EM«Ae 2 mN "em^n 1 4 x 10~ 41 cm 2 (50 GeV) 4 
a ~ 647T -. = — — — ~ — — 4 (5.23) 

where m^ is the mass of the nucleon. In order to fix a ~ 10~ 40 cm 2 , one must therefore 
choose < eDy/gA- For natural values of the Yukawa couplings, m^M < ^a, and so 
from Eq. 5.6 and the requirement that ~ 100 GeV follows the constraint that 

— > < 100GeV > 2 . (5.24) 



9a ~ D 



Y 
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Bounds on kinetic mixing [24] impose further constraints, requiring e < 0.005 for ~ 
1 GeV. Finally, the scale of supersymmetry breaking is proportional to e 2 g\ and to get 
splittings of order 100 keV requires that e ~ g A ~ 0.005. 

Trying to match the CoGeNT/DAMA anomaly [23] with light inelastic dark matter 
is challenging in this specific model because it is difficult to generate an 0(1O~ 38 cm 2 ) 
cross section. The primary tension arises from the mediation through the massive axial 
current. One approach could be to kinetically mix the vector current rather than the 
axial current. This does not suppress elastic scattering; however, in models of inelastic 
dark matter where the signal arises from downscattering, rather than upscattering, the 
elastic rate does not need to be suppressed. This possibility is not pursued further here 
but illustrates the rich phenomenology possible in composite dark matter models [20]. 

Following the above logic we choose the following MSSM parameters for all three 
benchmark models 





Mi 


Mec 


Auv 


50 GeV 


100 GeV 


900 GeV 


100 TeV 



(5.25) 



The parameters of the dark sector for the three benchmark points are chosen to be 
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Mixed 
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0.004 


30 GeV 
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1.4 


Heavy Scalars 


0.005 


0.5 


0.004 


30 GeV 
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1.0 
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The first two choices for g v cause £/(l)v to hit a Landau pole before the GUT scale. 
The Landau pole can be avoided by embedding ?7(l)v into a non-Abelian group, e.g. 
U(l) y C SU(2)v or C/(l)v x U(1) A C SO(4). The rather large values of e chosen 
here make for some tension with constraints from BaBar. These constraints may not 
apply to this model because A M may cascade decay through the Higgs sector before 
decaying into Standard Model particles [11]. These parameters lead to the following 
supersymmetric bound state mass scales 
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Supersymmetry breaking effects are encapsulated in the soft parameters 
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(5.28) 
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Figure 2: The ground state spectrum of the three benchmark models. In all cases the small, 
unlabeled splittings are of order B ~ 6 eV. 

The gauge-mediated contribution to rriy listed here is subdominant to gravity-mediated 
contributions, which give 

m -^%^leV (5.29) 
v M pl 



where yF susy ~ 100 TeV for these benchmark models. 

The resulting spectra are shown in Figure 2. The first parameter point, "Unmixed," 
realizes the scenario where B <C m so f t <C rnps- The various bound states have a 
mass 56 GeV and the lowest state o;_ is primarily scalar-scalar. The lowest states 
accessible by axial photon exchange are the pair of nearly degenerate scalars, uji + 
and w 2+ , which are heavier by an amount 5 = 64keV. For the second parameter 
point, "Mixed," the hierarchy of scales is instead B <C m so f t ~ mps and there is large 
mixing between the vector multiplet and hypermultiplet. The bound states have a 
mass 67 GeV and an iDM-compatible spectrum is again realized with several states 
available for upscattering. w_ is again the lightest state but the vector is now 
kinematically accessible. Because scattering to is velocity suppressed, however, the 
relevant splitting for iDM is 5 = 80keV between U- and wi + /w 2+ . For the third 
parameter point, "Heavy Scalars," the hierarchy of scales is instead B <C mps <C m so f t 
and mDM = 51 GeV. Because the hypermultiplet has large selectron components, the 
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hypermultiplet states move upwards and the spectrum is inverted, with (which is 
primarily fermion-fermion) and as the lightest states. The relevant splitting for iDM 
is 6 = 126 keV. 

From the results in Sec. 4.1 one can estimate the lifetime of the unstable states in 
the "Unmixed" benchmark spectrum. One finds ry ~ 10~ 18 sec for the vector multiplet 
and t% ~ 10~ 13 sec for the states in the hypermultiplet. Although the decay formulae 
from Sec. 4.1 are not directly applicable to the "Mixed" benchmark, the decay rates 
will be similarly fast, since ay is 0(1). 



6. Discussion 

The benchmark models discussed in Sec. 5.4 demonstrate that the non-relativistic su- 
persymmetric bound states discussed in this article can realize a wide range of bound 
state spectra with a rich hierarchy of scales. In any application of these models to 
dark matter phenomenology, a number of important issues must be addressed. In par- 
ticular one needs to examine constraints from direct and indirect detection as well as 
the implications for early universe cosmology. Although addressing these topics in any 
detail lies outside of the scope of this paper, in this section we briefly discuss some of 
the relevant physics. In the following our discussion is limited to the concrete model 
presented in Sec. 5. 

6.1 BBN Constraints 

In a realistic model where the dark matter sector is nearly supersymmetric, the hidden 
?7(l)v photon and photino that create the bound states will usually be relativistic at 
the temperatures relevant for BBN, T « 1 MeV. This means that they contribute 
to the energy density of the universe and will modify the successful predictions of 
standard BBN. Fortunately, the resulting perturbation is sufficiently small to be within 
observational constraints. This arises because the hidden sector is weakly coupled 
to the visible sector, typically e < 1CT 2 , and therefore the dark sector kinetically 
decouples [25] from the visible sector at a temperature above the GeV scale. The 
early kinetic decoupling reduces the number of effective degrees of freedom in the 
dark sector because many Standard Model degrees of freedom become non-relativistic 
between temperatures of 10 GeV and the QCD phase transition. This generally makes 
weakly coupled dark sectors with long range Coulombic interactions [26] safe from BBN 
constraints [27]. 
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6.2 Recombination 

If dark atoms are relevant for cosmology, then a significant fraction of the supersym- 
metric dark electrons and protons must recombine into supersymmetric atoms. The 
recombination of non-supersymmetric hydrogen-like dark matter atoms is studied in [9] . 
Supersymmetry adds new levels of complexity to the problem. First, supersymmetry 
introduces new processes where gauginos are emitted in de-excitation processes. Sec. 4 
discusses some of these de-excitation processes and, as a rule of thumb, gaugino emission 
processes are faster if the corresponding gauge boson emission process is a magnetic, 
spin-flip transition. For non-supersymmetric dark atoms, the most important processes 
for recombination are the Ly-a decay, the 2s double photon decay and the scattering 
process e~ +p + <r» H+j [9]. None of these are due to spin-flips and so gaugino emission 
processes are expected to be subdominant. 

Even if gaugino emission is subleading as expected, there are new electric transi- 
tions between the different, non-degenerate superspin levels with n = 2 and n = 1. This 
additional complexity makes the out of equilibrium problem harder to solve systemati- 
cally. The numerous new Ly-a lines corresponding to different superspin transitions are 
sufficiently degenerate that the Doppler broadening smears the energies of the states. 
For this reason photons from different transitions are indistinguishable implying that 
the optical depth of Ly-a photons is about the same as in the non-supersymmetric case. 
Although a detailed analysis is necessary, it appears that the physics of recombination 
in the supersymmetric case is parametrically the same as in the non-supersymmetric 
case. In particular for sufficiently large ay recombination is expected to be an efficient 
process. 

6.3 Molecules 

If supersymmetric atoms can form, it is possible that these atoms may further aggre- 
gate into supersymmetric molecules. This section briefly explores this possibility by 
examining the role that Bose/Fermi statistics plays in atoms and molecules. 

The ground state of regular diatomic hydrogen is a s = state, i.e. the spin wave- 
function for the electrons is antisymmetric under particle exchange and the spatial 
wavefunction is symmetric as it would be for a two-selectron atom. The statistics of 
the electrons thus does not affect the ground state spatial wavefunction and energy to 
lowest order. Supersymmetric hydrogen should therefore form diatomic molecules, and 
the binding energy is approximately the same as for non-supersymmetric hydrogenic 
systems [28]. 

In the Standard Model, further aggregation into molecules larger than H2 is pre- 
vented by the Pauli exclusion principle, which forbids more than two electrons from 
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being in the same orbital. In supersymmetric atoms, this aggregation is not forbid- 
den by Pauli because electrons can convert into their scalar superpartners, selectrons. 
Supersymmetric bound states will share orbitals more effectively and hence are bound 
more strongly. For non-relativistic molecules composed of N bosonic constituents, the 
binding energy scales as 



in contrast to molecules with fermionic constituents, whose binding energy grows lin- 
early with the number of constituents [29-31]. This means that macroscopic bound 
states formed from scalar constituents have an enormous binding energy. This suggests 
a scenario where a fraction of the dark atoms condense into huge "molecules," which 
may impose additional constraints on supersymmetric atomic dark matter if the for- 
mation of macro- molecules is too efficient. One possibility is that the supersymmetric 
macro-molecules could drive formation of microscopic black holes. 

6.4 Dark Matter Genesis 

In order for dark atoms to be a sizeable fraction of the Universe's dark matter, there 
needs to be a chemical potential for U(l) e+P generated in the early Universe [32]. The 
most compelling mechanism for generating an asymmetry in the dark matter sectors 
links the dark matter number density to the baryonic or leptonic number density, 
for recent work see references in [33]. The primary novelty with generating a chemical 
potential for the composite sector is that the minimal gauge invariant operator sourcing 
U(l) e+P is a dimension 2 operator, e.g. EP. This typically means that constructing 
a renormalizable Lagrangian requires more structure than for elementary dark matter. 
Normal inelastic dark matter requires ni>/n e+p > 1 and this means that if the e + p 
asymmetry is to be directly linked to the baryon asymmetry, then the interactions that 
equilibrate the chemical potentials in the dark sector and Standard Model must freeze 
out when the dark matter is relativistic. Alternatively, the e+p asymmetry should be 
generated through a mechanism like the out-of-equilibrium decay of a heavy particle. 

7. Conclusions 

Composite dark matter offers a rich phenomenology that has only barely been ex- 
plored in comparison to models in which dark matter is an elementary particle. Non- 
relativistic bound states offer one general class of composite dark matter models, and 
these typically involve a new mass scale that is incorporated either by hand or through 
a Higgs mechanism. In the later case the Higgs mass is radiatively unstable and su- 
persymmetry is a natural way to stabilize its mass. Supersymmetry breaking can be 
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weakly communicated to this sector, particularly if the interactions of the dark sector 
with the Standard Model are the dominant link to the supersymmetry breaking sector. 
If this is the case, then the composite dark matter will form nearly supersymmetric 
multiplets and the phenomenology of these states can be radically different from the 
non-supersymmetric case. 

Atomic Inelastic Dark Matter [9] was proposed as a model of inelastic dark mat- 
ter where the hyperfine splitting of the ground state of a hydrogen-like sector is the 
origin of the inelastic mass splitting. This article created a supersymmetric version 
of Atomic Inelastic Dark Matter that has many different features from the original, 
non-supersymmetric model. Most notably, the hyperfine splitting disappears in the 
supersymmetric limit and the ground state typically contains a scalar-scalar compo- 
nent. This introduces new interactions between dark matter and the Standard Model. 
The model accommodates spectra and scattering cross sections compatible with iDM 
phenomenology. 

This article has also constructed tools to help in the study of quasi-perturbative 
supersymmetric bound states and in incorporating supersymmetry breaking into the 
bound states. These tools illustrate that the form of the weakly broken spectrum and 
the composition of the various states is in many cases dictated entirely by supersym- 
metry at leading order. Supersymmetry also imposes strict restrictions on the allowed 
interactions, which simplifies the matching of effective interaction operators. More gen- 
erally supersymmetric bound states offer a rich laboratory for studying supersymmetric 
dynamics and interactions. 
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